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Dynamic equations are presented for polar binary mixtures containing ions in the 
presence of the preferential solvation. In one-phase states, we calculate the dynamic 
structure factor of the composition accounting for the ion motions. Microphase sep- 
aration can take place for sufficiently large solvation asymmetry of the cations and 
the anions. We show two-dimensional simulation results of the mesophase formation 
with an antagonistic salt, where the cations are hydrophilic and the anions are hy- 
drophobic. The structure factor S(q) in the resultant mesophase has a sharp peak 
at an intermediate wave number on the order of the Debye-Hiickel wave number. As 
the quench depth is increased, the surface tension nearly vanishes in mesophases due 
to an electric double layer. 



I. INTRODUCTION 



Much attention has been paid to the phase transition behavior arising from the Coulomb 
interaction among charged particles in various soft matters including electrolytes, poly- 
electrolytes, and gels [HO [3]. However, in most of the theoretical literature, the ion-dipole 
interaction has not been explicitly considered, which gives rise to a complex structure around 
each ion, called the solvation (hydration) shell, composed of several solvent molecules (those 
of the more polar component in a mixture solvent) jlj. The resultant solvation chemical po- 
tential fj, J sol depends on the ion species j and typically much exceeds the thermal energy T. It 
should also strongly depend on the composition for binary mixtures and the polymer volume 
fraction for polymer solutions, so it cannot be neglected at phase transitions and around 
composition heterogeneities. Recently, including the preferential solvation effect, several 
theoretical groups have begun to investigate the ion effects in electrolytes [HI El El El E], 
polyelctrolytespj)], and ionic surfactants [TT] . 
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We mention some experiments directly related to our theory. First, many authors have 
long observed salt-induced phase separation or homogeneization in aqueous binary mixtures, 
where the phase behavior is strongly altered even by a small amount of an salt |12j . Sec- 
ond, we mention a number of observations of salt-induced aggregates in near-critical binary 
mixtures [131 EH [IS], where the cations and anions are both hydrophilic. In one-phase 
states [131 [Hj, heterogeneities extending over a few micrometers have been detected by light 
scattering with addition of a salt (for example, 17 mass% NaBr in mixtures of H 2 0+3- 
methylpyridine (3MP) [14]). In two-phase state [IS], a macroscopic thin plate has been 
observed at a liquid-liquid interface, which presumably consists of aggregates of solvated 
ions. Third, we mention recent small-angle neutron scattering experiments by Sadakane et 
al. [El |T7j. They added sodium tetraphenylborate NaBPli4 at 100 mM to a mixture of D2O 
and 3MP to find a peak at an intermediate wave number g m (~ 0.1 A -1 ). The peak height of 
the SANS intensity was much enhanced with formation of periodic structures. Their salt is 
composed of hydrophilic Na + and hydrophobic BPhJ. Furthermore, the mixture exhibited 
colors changing dramatically on approaching the criticality at low salt contents (~ 10 mM). 

Hydrophilic and hydrophobic ions interact differently with the composition fluctuations 
in mixtures of water+less polar component. They behave antagonistically in the presence 
of the composition fluctuations. We may predict formation of a large electric double layer 
at liquid-liquid interfaces much reducing the surface tension and formation of mesophases 
for sufficiently large solvation asymmetry [SI E]. However, we do not know the details of 
the phase transition of binary mixtures with an antagonistic salt. Originally, Nabutovskii 
et al. [IE] pointed out a possibility of mesophases in electrolytes assming a coupling between 
the composition and the charge density in the free energy. 

In Section 2, we will present a short summary of the statics of binary mixtures containing 
ions accounting for the preferential solvation. In Section 3, dynamic equations for such 
systems will be given and, as an application, the dynamic scattering amplitude will be 
calculated. In Section 4, we will then numerically examine the mesophase formation induced 
by antagonistic ion pairs. 
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II. GINZBURG-LANDAU FREE ENERGY 
A. Electrostatic and solvation interactions 

We consider a polar binary mixture containing a small amount of salt. The composition of 
a water-like component is written as <fi. The cation and anion densities are written as n\ and 
i%2 with charges Z\e and Z 2 e. In the monovalent case we have Z\ — 1 and Z2 = — 1. They 
are sufficiently dilute and their volume fraction is negligible. The chrage density is given 
by e(Z\n\ + Z 2 n 2 ). The variables 0, n\, and n 2 vary smoothly in space. The Boltzmann 
constant will be set equal to unity. As the geometry of our system, our fluid is between 
parallel metallic plates in the region < z < L. The lateral dimensions in the xy plane are 
much larger than L. The surface charges on the plates can give rise to an applied electric 
field. In the following theory, we fix the charges on the plates such that their electrostatic 
energy is kept fixed [T9] . 

The Ginzburg-Landau free energy functional of our system is written as F = J drf with 
the free energy density [5J E] , 

/ = /o(0,T) + ^|V0| 2 + £f;i 
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kifajV 0) - 1 - gj(f> 



(2.1) 



The first two terms constitute the usual Ginzburg-Landau free energy density. The chemical 
part fo = fo((p,T) depends on and T and the coefficient C of the gradient term is of 
the order a 2 ~ d in d dimensions, where a is the molecular radius. The third term is the 
electrostatic free energy, where E = — V$ is the electric field and the electrostatic potential 
$ satisfies the Poisson equation 

- V-£(0)V$ = 47re(Zini + Z 2 n 2 ). (2.2) 

The dielectric constant e((p) can depend on the composition 0. In our previous work the 
linear composition dependence 

e{4>) = e + ei<p (2.3) 

has been assumed, where e$ is the dielectric constant of the less polar component and 
£o + £i is that of the water-like component. In such cases Eq > and E\ > 0. This linear form 
approximately holds in some polar binary mixtures [20]. The last term in (2.1) consists of the 
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entropic part and the solvation contribution of the ions, where the parameters gj represent 
the solvation strength. The choice of the volume vq is arbitrary and is taken to be the solvent 
molecular volume (see (5.1) below). In this work we neglect the image interaction arising 
from inhomogeneous dielectric constant or from nonvanishing E\ in our theory (5j EJ |2T] . 
The interfacial ion distribution is then produced by the preferential solvation among the 
ions and the mixture solvent. The image interaction is weakened with increasing the salt 
density and/or approaching the critical point. 

With (2.1) we may calculate the chemical potentials h = 8F/S(f) and fij = SF/Srij. They 
are written as 

f = f -^VV-f^-E^, (2.4) 

j 

^ = ln(n^o) - 9j4> + yZje®, (2.5) 

where f' Q = dfo(4>)/d(j) and E\ = de/dcp. If the system is in equilibrium, h, /ii, and ^ are 
homogeneous constants. When the system undergoes a macroscopic phase separation with 
a planar interface separating polar and less polar regions, we may calculate the interface 
profiles of the composition and the ions [6]. In equilibrium the composition difference A0 
and the potential difference A$ satisfy 

eA$ = T{g 1 -g 2 )A<f>/{Z 1 + \Z 2 \), (2.6) 

from the charge neutrality in the bulk regions. The A$ is called the Galvani potential 
difference in electrochemistry [221 123] . 

The solvation free energy may be written as / so i = J2j l J lo\( < P) n ji where fi J so \((fi) is the 
solvation chemical potential of the ion species j. It is assumed to depend on as 

/C(0) = /4oi(O)-^0. (2.7) 

Here the first term in the right hand side gives a contribution linear in nj in / so i and is 
not written in / in (2.1), while the second term yields the solvation coupling terms in / 
between the ions and the composition. We remark on the magnitude of gj. In aqueous 
mixtures, it is positive for hydrophilic ions and negative for hydrophobic ions. In two-phase 
coexistence, the difference of the solvation chemical potential between the two phases is given 
by A/ig ol = TgjAcf), which is identical to the standard Gibbs transfer free energy (per particle) 
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in electrochemistry J22J [23]. Data of Afi 3 sol are available for water-nitrobenzene at room 
temperatures in strong segregation (where A<p = 1). For example, A/jp^/T = gjAcf) = 13.6 
for Na + , 15.3 for Li + , 26.9 for Ca 2+ , 11.3 for Br~, and 7.46 for I - as examples of hydrophilic 
ions, while it is —14.4 for BPhJ (tetraphenylborate) as an example of hydrophobic ions. 
The anion BPhJ consists of four phenyl rings bonded to an ionized boron, acquiring strong 
hydrophobicity. Note that Sadakane et al. [161 E] used NaBPli4. Thus the preferential 
solvation effect can be very strong. However, it has mostly been neglected in theories of 
electrolytes and soft matters, though it strongly influences phase transitions in such systems. 

When phase-separation occurs macroscopically, a liquid-liquid interface appears. If the 
space dependence is along the z axis, the surface tension is expressed as \±0\ 

7 = 2 j dz(f g -f el ), (2.8) 

where f g = C\ V0| 2 /2 is the gradient free energy density and f e \ = eE 2 /8n is the electrostatic 
free energy density. Up to linear order in the ion densities, we may also derive the expression, 

7 7o - TT + leh (2.9) 

where 70 is the surface tension without ions, Y is the surface adsorption of ions, and 7^ is 
the electrostatic contribution given by 7 e ; = — / dzfei < 0. For antagonistic salts with large 
\gi\, \j e i\ is much amplified due to the electric double layer at the interface. 



B. Structure factor in one-phase states 

In our previous papers [3 |6], we examined the structure factor S(q) = (\4>q\ 2 ) of the 
composition fluctuations with wave number q = \q\ in one-phase states with salt, where 4>q 
is the Fourier component of the composition deviation 5<j)(r) = <j)(r) — (0) with wave vector 
q. Hereafter (■ ■ ■) denotes the thermal average. We readily obtain S(q) if the fluctuation 
contributions to F are calculated in the bilinear order. The resultant free energy part 
is written as 5F. The thermal fluctuations obey the Gaussian distribution oc e~ SF ^ T in 
equilibrium in the mean-field theory. Hereafter we consider the monovalent case Z\ = 
—Z 2 = 1, where the average ion densities are written as (ni) = (n 2 ) = n e . 

From (2.1) some calculations give 



T 
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i 1 2 i* 



2n 



(2.10) 



where n^g and /?g are the Fourier components of rij(r) and p{r) = nx(r*) — ^(t*), respectively, 
and = e 2 /eT is the Bjerrum length. We define 

f = 9 2 /o(0)/^ 2 . (2.11) 

The average composition (0) is simply written as 0. Here we may treat e as a constant when 
we treat the small thermal fluctuations. By minimizing 5F with respect to njq at fixed 4>q, 
we obtain 5F/T = Y,q \<Pq\ 2 /2S{q) with 



2,„2 
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1 ^— - 

k 2 + g 2 



(2.12) 

where k = (SvrfB^) 1 / 2 is the Debye wave number and the parameter 

lp = (lQnC£ B )- 1/2 \ gi -g 2 \ (2.13) 

represents asymmetry of the solvation of the two ion species. The structure factor thus 
obtained is analogous to that for weakly charged polyelectrolytes pJJ [2U [25] . 

The second term in the right hand side of (2.12) gives rise to a shift of the spinodal curve 
[T2] . For example, if the cations and anions are hydrophilic and g\ ~ ~ 15, the shift term 
is of order — 500n e and its magnitude can be appreciable even for v^n e <ti 1. On the other 
hand, 7 P can be increased for antagonistic salts composed of hydrophilic and hydrophobic 
ions [3 El [TH [T7]. From the last term in (2.12) a Lifshitz point appears at j p — 1. For 
7 P > 1, S(q) exhibits a peak at an intermediate wave number q m . Since the derivative of 
the right hand side of (2.12) with respect to q 2 vanishes at q = q m , we find 

Qm = (7 P - 1) 1/2 ^- (2.14) 

The peak height is given by S(q m ) = l/(f — r m ), where 

r m = (^i+^) 2 y + C(7p-l)V. (2.15) 

For f < r m , mesophase formation takes place, as will be studied in Section 4. 
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III. DYNAMICS 

A. Dynamic equations for composition, ions, and velocity 

We present the dynamic equations for <fi, rii, ri2, and the velocity field v [26]. The fluid 
is assumed to be incompressible and isothermal. That is, we require 

V-v = (3.1) 

and treat the mass density po and the temperature T as constants. Then <ft and rij obey 

^ + V-(<H = LoV 2 ^, (3.2) 
9Uj + V ■ (rijv) = DjV ■ rijV^ 



(3.3) 



Dp. 



Zje 

Vrij - -T^rijE - gjUjVcp 



where h and and \ij are given in (2.4) and (2.5), L is the kinetic coefficient (with L /vo 
representing a diffusion constant), and Di and D 2 are the ion diffusion constants in the 
solvent. The momentum equation is expressed as 

Po m = ~ Vpi ~ V-n + r/ V 2 v, (3.4) 

The first term on the right hand side ensures the incompressibility condition (3.1) and p\ 
satisfies 

V 2 pi = -E V « V ^> (3-5) 

a/3 

where V a = d/dx a with x a = x,y, z. We introduce the reversible stress tensor II = {n Q( g} 
(a, (3 = x, y, z) in the form, 

n a/3 = TCV a <pV p <P - ^-E a Ep. (3.6) 

47T 

where the first term is well-known in critical dynamics [26] and the second term is a part of 
the Maxwell stress tensor (with its diagonal part being included in pi) [27J. 
We determine II from the relation, 

V ■ n = <f)Vh + njVjJLj. (3.7) 
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If the above relation holds, the total free energy Ft = J dr[f + poV 2 /2] including the fluid 
kinetic energy changes in time as 

d r 

—F T = - J dr[e <j> + e vis + e ion ], (3.8) 
where the terms in the brackets are the heat production rates in the bulk given by 

= L Q \ Vh\ 2 , e vis = VoYl I v ^/?| 2 , 

a/3 

eion = E^|V/i,r/^ (3-9) 

3 

The surface terms are omitted in (3.8). Owing to dF T /dt < 0, the system tends to equilib- 
rium if there is no externally applied flow. 

In our dynamic equations we neglect the random source terms [26J, which are related 
to the transport coefficients Lq, Dj, and i]q via the fluctuation-dissipation relations. They 
are needed to describe the dynamics of the thermal fluctuations and to calculate the time 
correlation functions. 



B. Stokes approximation 

Without macroscopic flow, the viscous motion of v is much faster than the diffusive 
motions of <ft and rij. Here Lq/vq, Di, and D 2 are estimated by the Stokes formula (Dj ~ 
T/67rr] aj with aj being the molecular size), so they are much smaller than the kinematic 
viscosity r^o/po- Then we may well neglect the acceleration of the velocity in (3.4) to obtain 
[26] 

v a (r) = fdr'J2 T a p{r - r')X p {r') (3.10) 

p 

where X a (r) = —^ZpVpU. a p{r) is the force density acting on the fluid and T a p{r) is the 
Oseen tensor. This Stokes approximation has been used in numerical analysis of spinodal 
decomposition in the literature j2H]- The free energy F = J drf changes in time as dF/dt = 
— J drle^ + e vis + e ion ] < as in (3.8), where e vis is replaced by 

a 

Here / dre v - ls > from the expression (3.10). 
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C. Ionic local equilibrium 

The composition evolution can be much slower than the ionic motions particularly near 
the critical point. In such cases, the ion distributions are expressed in terms of <f> and $ as 

rij = n] exp{g j( f) - ZjeQ/T), (3.12) 

where the coefficient n°- is determined from the conservation of the ions / drrij(r,t) =const. 
In numerical analysis this approximation is convenient to examine the mesophase formation 
for large gj. 



IV. RELAXATION OF THE THERMAL COMPOSITION FLUCTUATIONS 



A. Time-correlation function 



In this section, we calculate the time-correlation function of the Fourier components of 
the composition fluctuations, 



G{q,t) = (<f> q (t)M0y), 



(4.1) 



in one phase states. This function can be measured by dynamic scattering. It is of interest 
how it relaxes on approaching the spinodal point and how it is influenced by the ion diffu- 
sion. The thermal hydrodynamic fluctuations are governed by the linearized hydrodynamic 
equations of (3.2) and (3.3) with random source terms added. That is, they obey linear 
Langevin equations [26]. In this section, without explicit introduction of the noise terms, 
we will calculate the time-correlation functions of the form (Aq(t)(f)q(0)*) with t > 0, where 
A — 4>, Hi, an d n 2 . We also assume that the cations and the anions have the same diffusion 
constant or D 2 — Di, which much simplifies the calculation. 
From (3.3) G(q,t) obeys the linear equation, 





Of 



+ r (q) 



G = L q 2 



giGi + #2^2 



(4.2) 



where G(q,t) is written as G and T (q) is the decay rate without ions, 



T (q) = L q 2 (f + Cq 2 ). 



(4.3) 
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Here we write G\ = (niq(t)4> q (0)*) and G 2 = (n 2 q(t)4>q(0)*}, where n\q and n 2 q are the 
Fourier components of n\ and n 2 . The equations for the combinations G\ ± G 2 read 

d 



dt 



d_ 
dt 



+ D 1 (q 2 + K 2 ) 



(G 1 + G 2 )=n e (g 1 + g 2 )D 1 q 2 G, 
(G 1 -G 2 )=n e (g 1 -g 2 )D 1 q 2 G. 



(4.4) 
(4.5) 



Use has been made of the fact that the Fourier component of the electric potential is $g = 
4-7re(niqr — n 2 q)/eq 2 from the Poisson equation (2.2), where the dielectric constant e may 
be treated as a constant. The convective terms in (3.2) and (3,3) vanish in the linear order 
without velocity gradient. 

It is convenient to calculate the Laplace transformation 



J 



dte- nt G(q,t). 



(4.6) 



The Fourier transformation I(q, u) = /f^ dte tujt G(q, t) is related to G(q, £1) by 

I(q,u) = 2Re[G(q,icu)], 



(4.7) 



where Re[- • •] denotes taking the real part. Some calculations give the following expression, 



G(q,Q) = S(q) 



L q 2 /S(q) 



The ionic correction Z(q, Q)(oc n e ) depends on q and Q as 



(4.8) 



Z(q,n) = ^n e L q 2 



(9i + 92f 



n + D iq 2 

(9i - 92?q 2 



[n + D 1 (q 2 + K 2 )}(q 2 + K 2 )_ 
In deriving (4.8) and (4.9) use has also been made of the static relations, 

((n iq + n 2 q)4>* q ) = n e (g x + g 2 )S(q), 

(1 + K, 2 /q 2 )((n iq - n 2q )(p q ) = n e (g 1 - g 2 )S(q), 



(4.9) 



(4.10) 



which follow from (2.10). These equal-time correlation functions appear in the Laplace 
transformation of S(q, t) in the presence of the random source terms. The presence of 
Z(q,Q) in (4.8) makes the relaxation of G(q,t) complicated. 
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B. Relaxation near the spinodal point 

We obtain the exponential relaxation, 

G{q,t) S{q)e~ r ^ t , (4.11) 

near the spinodal point. Here the decay rate T(q) is assumed to be much smaller than -Dig 2 . 
Then we may set Q = in Z(q, Q) to find 

m = Loq2 / s % (4.i2) 

W l + n e B( q y K } 



where B(q) = Z(q,0)/n e is written as 

B(,) = I 



(gi + 92? {gi - ftQV 1 , , 

2 2(g 2 + K 2 ) 2 J' 1 J 

If \gi\ and l^l ar e very large, the ionic correction n e B(q) can be noticeable even for t> n e "C 1. 
For 7 p > 1, T(g) tends to zero first at q = q m on approaching the spinodal point. 

C. Long wavelength limit 

In dynamic light scattering experiments, we should consider the long wavelength limit, 
where we set q < k, T (q) = D^q 2 , and Z(q, SI) = aDiq 2 /(Q + D x q 2 ) with 

D<j> = L f, (4.14) 
a = L ( 9l + g 2 ) 2 n e /2D 1 . (4.15) 

Here = lim^o F(q)/q 2 is the diffusion constant of the composition in the long wavelength 
limit without ions. The dimensionless parameter a is proportional to n e and increases steeply 
with increasing n e for g\ + g 2 ^> 1. In dynamic light scattering without ions, tends to 
zero near the critical point (being given by the Kawasaki formula T/6in] C, with £ being the 
correlation length) [26] . In this limit we obtain 

G(q,0) tt + {l + a)D iq 2 

S(q) ~ tt 2 + Vq 2 tt + D X D^ 

" + (4.16) 



tt + D_q 2 Vt + D+q 2 

where T> — (1 + a)Di + in the first line. The two new diffusion constants D_ and D + in 
the second line are expressed as 

V 



D ± = -±-jV*-4D+D 1 . (4.17) 
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The partition coefficient (3 is of the form 

^-2 + 2(D + -D_) • (418) 
The inverse Laplace transformation of the second line of (4.15) yields the time-correlation 
function exhibiting a double-exponential decay, 

°Ml = (3e- D -« 2t + (1 - (3)e- D ^ H . (4.19) 
S(q) 

(i) For very small ion concentrations there can be the situation where a C 1 and D\a <C 
D+<D X . In this case we have ZL = £> and D + ^ Dj with (3=1- aD 4> j{D x - D^) 2 + ■■■. 

(ii) We may suppose the case D s <C D\a . In this case we have D- = D^/(l + a), D + = 
£>i(l + a), with (3=1- [(2a - 1)/2(1 + a) 2 ]ZV£i + • • •. 

V. SIMULATIONS AT THE CRITICAL COMPOSITION 

We numerically examine phase ordering with a strongly antagonistic salt at the critical 
composition (0) = 1/2. The spatial dimensionality d has been equal to three so far. However, 
we here present preliminary simulation results in two dimensions. 

A. Numerical method 

In our simulation, we choose the chemical free energy density fo in (2.1) in the Bragg- 
Williams form, 

^/o = 01n0 + (1 - 0) ln(l - 0) + X 0(l - 0), (5.1) 

where v = a d is the solvent molecular volume and x is the interaction parameter dependent 
on T. The parameter f in (2.9) is given by f = [1/0(1 — 0) — 2x\/v . Space and time will 
be measured in units of a and 

to = v a 2 /L , (5.2) 

where L Q is the kinetic coefficient in (3.2). Integration of the dynamic equations is performed 
on a 256 x 256 square lattice, so the system is in the region < x, y < 256a. Supposing the 
monovalent case, we set 

v C = a 2 , g x = -g 2 = 15, 

n n o _ 0.16a 4 
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Then we obtain 7 P = 2.44 from (2.11) and mesophases are realized with increasing \. These 
values of g\ and gi are realistic in view of the data of the Gibbs transfer free energy, as 
discussed below (2.7). The correlation length is defined by £ = [Cv /\A — 2x|] 1//2 , which is 
equal to a for \ — 2.5. 

The velocity field v is determined by the Stokes approximation (3.10) and the ion densities 
ri\ and n 2 by the Poisson-Boltzmann expressions (3.12) (the latter being justified in the limit 
Dj — > oo). In the dynamic equation (3.2) for cf) we put a random source term to calculate 
the structure factor, 

^ + V-(^) = L V 2 ^-V-^, (5.4) 
Here j R is the random diffusion flux characterized by 

(jRa(r,t)j Rt3 (r',t')) = 2L5 a/3 5(r - r')5(t - t'). (5.5) 

where a, (3 = x, y. The noise strength L should be equal to the kinetic coefficient Lq to 
ensure the equilibrium distribution (oc e~ F//T ). In this paper, however, we set L = 10~ 8 Lo to 
detect the composition patterns unambiguously. In one phase states, <fi remains nonvanishing 
due to j R , yielding a structure factor proportional to the mean field structure factor S(q) 
in (2.12), where the proportionality constant is L/L = 10 -8 (not shown here). In two 
phase states, j R serves to trigger phase ordering, yielding a structure factor composed of the 
domain contribution. The same structure factor follows even if we set j R = in the course 
of domain growth. It is worth noting that the random source terms are mostly neglected in 
the literature of phase ordering dynamics [26J . 

In our simulations we start with the initial condition 0(r,O) = 1/2 at t — 0. Small 
disturbances of <ft are subsequently produced by the small random flux j R in (5.4), which 
grow into patterns in two-phase states. For g\ = —g 2 and (0) = 0.5, use of (2.15) yields the 
linear instability criterion, 

2-X< l -Ca\ lp -l) 2 n\ (5.6) 

where the right hand side is 78fon e from (5.3). Hereafter the Debye wave number is k = 
8.7rig/ 2 with n e = (ni) = (n 2 ) being the average ion density. 
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B. Mesophase formation in shallow quenching 

Here we study the phase ordering at the solvent criticality x = 2 and ((f)) = 1/2, where 
instability occurs for n e > 0. In Fig.l, we show the time evolution of a normalized char- 
acteristic domain size 27r/aq p (t) for various n e . In terms of the time-dependent structure 
factor S(q,t) = (\<t>q(f)\ 2 ) we define 

<&»(*)=£<?£(?, (5-7) 
q q 

In Fig.l, q p (t) tends to a constant expressed as 9.22n^ 2 at long times. It nearly coincides 
with q m = (7 P — l) 1 / 2 /? = lOAn^ 2 in (2.14). In Fig.2, the steady-state structure factor S(q) 
is given for three ion densities, where all the curves arise from the domain structure and are 
not affected by the small noise term in (5.4). Our S(q) exhibits a sharp peak at q = q m and 
a second peak at q — 3q m . The peak height at q = q m is a constant of order unity nearly 
independent of n e . This can be explained as follows. It is known that a domain structure 
gives the structure factor of the form S(q) = (A<p) 2 £ d S*(q£), where A0 is the composition 
difference between the two phases, £ is the domain size, and S*(x) is a scaling function. In 
our two-dimensional case, we have A0 oc n l J 2 and £ ~ 2-k /q m oc n~ 1 / 2 , so S(q) ~ S*(q£). 

In Fig.3, we display <f>(r,t) and ni(r*,t) at t = 3000t for v Q n e = 0.0005 (left) and 0.0025 
(right). In Fig.4, we present cross sections of 0, n\, ri2 in the upper panel and those of 
the gradient free energy f g = TC\V(p\ 2 /2 and the electrostatic energy f e = sE 2 /8ir, and 
their difference in the lower panel. These quantities vary mildly without sharp interfaces as 
functions of x at fixed y = 64a. We notice that the difference f g — f e is small. In Fig. 5, 
their space averages, (f g ) = J drf g /V and (f e ) = J drf e /V, are demonstrated to be nearly 
the same at long times, where V is the system volume. 

We argue why (f g ) = (f e ) holds in steady states in weak segregation. If the ion density 
is small at shallow quenching, the composition is weakly segregated and is composed of the 
Fourier components with q = \q\ = q m . As in the weak segregation case of block copolymers 
[29] , the deviation 5<fi — <f) — (0) is expressed as 

6<j ) = Y,A q e^ r , (5.8) 

q 

where the coefficients Aq are sharply peaked at q = q m . With this form, the space average 
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t/L 



FIG. 1: Characteristic domain size 2ir / 'q m (t)a vs time at the solvent criticality for von e 
0.001,0.0015,0.002, and 0.003. The saturated value of q p (t) nearly coincides with q m in (2.14). 




FIG. 2: Steady-state structure factor S(q) of the composition for von e = 0.001,0.002, and 0.003, 
where the the solvent is at the criticality {x = 2 and {</)) = 1/2). 



of f g is written as 



T 



C 
2V 



q 



'■)■ 



(5.9) 



where (5<f) 2 ) = J2q \Aq\ 2 /V. Linearlizing (2.2) and (3.12) with respect to 5(f) in the monova- 
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n,v ft =0.0005 



»,v =0.0025 




0.005 



FIG. 3: Patterns of (f>(r,t) (top) and ni(r,t) (bottom) at t = 3000i for n e v = 0.0005 (left) and 
0.0025 (right) at the solvent criticality. 



lent case, we obtain the electric potential [10J, 



E 



2e ^7? q 2 + k. 2 



Age**, 



From (2.13) the average electrostatic energy is written as 

C/e) = 1 v C^KW 

T 2V ^ (q 2 + K 2 ) 21 q 



2(ql + K 



2\2 



'I- 



From q m = ( 7p - lf' 2 K in (2.14), we find (f g ) S (/ e ). 



(5.10) 



(5.11) 



C. Mesophase formation in deep quenching 



Next we examine the case of deep quenching by setting x — 2.5 with (0) = 0.5, where 
the interface thickness is £ = a. In Fig.6, we show the time evolution of the characteristic 
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FIG. 4: Cross sections of cj>, VqUi and von2 (top) and those of f g vo/T, —f e vo/T, and (f g — f e )vo/T 
(bottom) for n e vo = 0.0025 in the region 100 < x < 180 at y = 64, where the solvent is at the 
criticality. Use is made of the data producing the right images in Fig. 3. 

domain size 2ir/q p (t), where q p (t) is denned by (5.7). For v§n e = 0.0025 the domain size 
approaches a constant, while for t> n e = 0.0015 and 0.0005 its growth still continues in the 
simulation but is extremely slow at the end of the simulation (t/to = 6000). In Fig. 7, the 
structure factor S(q) is shown for v§n e = 0,0.0005,0.0015, and 0.0025. The structure factor 
around the peak is of order 100t>o and is much larger than the thermal level. 

In Fig.8, we display (j)(r,t) and ni(r*, t) at t = 3000t for v n e = 0.0005 (left) and 
0.0025 (right). As a marked feature for v n e = 0.0005, the cations (anions) are confined 
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FIG. 5: Space averages of f g v /T = v C\V(j)\ 2 /2 and f e v /T = v e\ V$| 2 /8vrT vs time t for 
n e vo = 0.0025, where the solvent is at the criticality. 



T 




| h„v„=0.0025 

°0 2000 4000 6000 

t/t 

FIG. 6: Characteristic domain size 2tt / 'q m (t)a vs time in deep quenching with \ = 2.5 and 
((f)) = 1/2 for v n e = 0,0.0005,0.0015, and 0.0025. 

in the water-rich (water-poor) regions. Because of the small ion density here, the ions 
change discontinuously at the interfaces and are homogeneously distributed in the preferred 
domains. On the other hand, for f n e = 0.0025, the ions are localized near the interfaces. In 
Fig.9, we show cross sections of 0, rii, n 2 (top) and those of f g , —f e , and f g — f e (bottom). We 
can see electric double layers at the interfaces in accord with the theory [6]. The difference 
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FIG. 7: Structure factor S(q) in deep quenching with \ = 2.5 and ((f)) = 1/2 for v$n e = 
0,0.0005,0.001,0.0015, and 0.0025. The numbers in the figure denote n e v . 

fg ~ fe turns out to be small in steady states. In Fig. 10, we demonstrate that their space 
averages nearly coincide at long times. From (2.8) we recognize that the surface tension 7 
nearly vanishes in steady states. 




VI. SUMMARY AND CONCLUDING REMARKS 



In this work, we have presented dynamic equations for binary mixtures containing ions, 
where the free energy includes the solvation interactions, (i) As the first application, we have 
calculated the dynamic structure factor G(q, t) in (4.1) in one-phase states accounting for the 
ion motions. Its relaxation is slowed down on approaching the spinodal as in (4.12), which 
occurs at an intermediate wave number q m for 7 P > 1. Here 7 P is the asymmetry parameter 
of solvation. It also exhibits a double-exponential relaxation in the long wavelength limit (in 
dynamic light scattering) as in (4.19). (ii) As the second application, we have numerically 
demonstrated emergence of mesophases with addition of an antagonistic salt, though our 
simulations are in two dimensions and at the critical composition. We have obtained a 
dramatic increase of the structure factor S(q) at an intermediate wave number in Figs. 2 
and 7 in accord with the experiment of Sadakane et al.[W\ [T7] . In these mesophases, the 
gradient free energy and long-range electrostatic energy are balanced as in Figs. 5 and 10. 
We have found that the surface tension vanishes in the mesophase in deep quenching. 
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n e v =0. 00001 « t ,v =0.0025 




o^m^m ■ 0.00003 om 10.0005 



FIG. 8: Patterns of <j)(r) (top), n\{r) (middle) and 712(f) (bottom) in deep quenching at t = 3000io 
for n e v = 0.00001 (left) and 0.0025 (right), where x = 2.5 and (4>) = 1/2. 

The present simulation is still very preliminary and more systematic analysis is needed 
in future work. In particular, the phase diagram in the parameter space of \i (0)) an d n e 
is required. While our simulation captures some salent features of the neutron scattering 
experiments [TO], [T7], the calculated structure factor S(q) cannot be compared with the 
observed intensity quantitatively. There are many parameters in our theory and we cannot 
judge whether or not our choice in (5.3) is appropriate for the experimental system. In 
particular, the solvation parameters gi and #2 are not known for mixtures of D2O and 3MP. 
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FIG. 9: Cross sections of <f>, VqUi and von2 (upper plate) and those of f g vo/T, —f e vo/T, and 
(f g — f e )vo/T (lower plate) for n e vo = 0.0025 in the region 100 < x < 180 at y = 64, where x = 2.5 
and ((f)) = 1/2. Use is made of the data producing the right images in Fig. 8. 

Our simulation suggests that addition of an antagonistic salt to a binary mixture can 
decrease the surface tension of a macroscopic liquid-liquid interface even to zero. We may 
then predict a salt-induced interface instability, leading to emulsification. We also men- 
tion measurements of the dynamic scattering, the electric conductivity, and the rheological 
properties, as new experiments using an antagonistic salt. 
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